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Heat flow and thermal profile in a ID harmonic lattice with coordinate-dependent masses has 
been calculated in the thermodynamic limit. It is shown in the particular example of a ID harmonic 
lattice with Unearly increasing masses that in standard Langevin conditions of contact, a temperature 
gradient can form, and Fourier's law can be obeyed. 



It is well known [1^ that in one-dimensional harmonic 
systems no temperature gradient can evolve, and neither 
can Fourier's law be obeyed, i.e., the dependence of heat 
flux J on the size of a system N cannot be described by 
the relation J ~ 1/A^. 

This conclusion was drawn on the basis of a consider- 
ation of only two systems for which these results can be 
obtained with analytical accuracy, more specifically, sys- 
tems of particles with equal or alternating masses [51, 15] . 
It came therefore as a surprise when numerical simula- 
tions demonstrated that in systems with linearly j4] or 
exponentially [S] increasing masses a temperature gradi- 
ent can form. 

Interest in ID systems in which mass depends on coor- 
dinates did not come as a surprise. It stemmed primar- 
ily from studies of functional-gradient materials (FGM) 
which can be employed in devices designed for monitor- 
ing and control of heat flow on a nanoscale |5H8j . 

The present study demonstrates a method of solution 
of linear equations describing such functional-gradient 
materials. In particular, for a system with an arbitrary 
dependence of masses on coordinate an exact analytical 
solution has been derived, which describes the tempera- 
ture proflle and the heat flow. It will be shown that in 
a system with a linear mass distribution heat gradient 
does indeed exist, and Fourier's law is indeed obeyed in 
certain cases. 

Consider a one-dimensional chain of + 1 particles 
which interact harmonically with their nearest neighbors 
with a spring constant K. The momentum and displace- 
ment n of a particle of mass m„ will be denoted by p„ 
and Xn, respectively. The Hamiltonian of this system 
(Fig. [T]) can be written in the form: 
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The equations of motion for such a system become: 

Xn COS /3n = ^ {Xn+1 + Xn-l) , 
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(2) 



where cos/3„ = (1 — w^m„/2if). 

In the thermodynamic limit the solution of this system 
of equations can be obtained in analytical form. 

To do this, we use Fourier's method by writing Xn in 
the form 
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FIG. 1. Schematic of a harmonic system with coordinate- 
dependent masses 



f{k) exp{ikn)dk. 
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We choose the limits of integration ki , k2 such that the 
function f{k) and all its derivatives tend to zero at these 
points. In this case we come to the equality 
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We further assume that the mass to„ and /3„ are functions 
of n/N. Now the function cos (3{n/N) can be expanded 
in the Taylor series with respect to the variable n/N, 
with Eqs. ^ converting into a differential equation for 
function /: 



/ = fcos{k), 
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Only in rare cases, such as a linear dependence of particle 
masses on index: 



Ml - Mo 

nin = Mq H — n, 
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this differential equation allows an exact solution. 

Because to this case we are going to revert more 
than once, we note that Eqs. p| are solved in terms 
of Hankel functions: x„ = H^^{z) , h':^^^(z) , where 
7 = zcos/3(0),z = N/{cosP{l) - cos/3(0)). This can 
be verifled by direct substitution using the recurrent 
relation 2n/zZ„(z) = Z„+i(z) + Z„_i(z), Z„(z) = 

In a general case, we are interested in the solution in 
the thermodynamic limit. We are going to look for the 
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solution to Eq. Q subject to the condition cxo in 

the form: 

fik) = ip{k)e's{k}N ^ 

In this case we come to the following relations for g and 
ip: 
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9 = - r'{z)dz, 



where /3~^ is a function inverse of /?. Note that (p does 
not depend on frequency. 

Next we choose the path of integration in Eq. ^ such 
that it will descend most steeply from the saddle point. 
Substituting the expressions for g and (p, we come in the 
limit that ^ cx) to a particular solution to Eqs. (p|: 



conditions of contact of the system under consideration 
with heat reservoirs. For the standard (Langevin) con- 
tact. SlC'^) — i^jSifiSjfi, S/j(a;) = iLO^Si^N^j.N ■ After 
a few straightforward transformations similar to those 
made in Ref. we come to: 
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\Do.N ~ ^{io){Di,N + Do,N-i) + S(u;)2i^i,w_ir 

(11) 

where Di „i is defined to be the determinant of the sub- 
matrix of K ~ uj'^M beginning with the Z-th row and 
column and ending with the m-th row and column. This 
determinant can be readily derived, because we know the 
general solution for Eqs. ([2|: 



(12) 



The solution of equation Dqn ~ with respect to fre- 
quencies uj yields the dependence of the wave vector k on 
the frequency of vibrations in the system under consid- 
eration k{u!): 



Vsin/3(a) 



4>{a)e 



ibl3(a) 



iN f f3{x)dx 



(8) 



where n is taken in the form n = aN + b, 4i{k) — 

We finally come to the general solution of Eqs. ^ in 
the form Xn = Ayn + -By^- 

Using the standard non-equilibrium Green function 
method [TJ [101 E] ^nd assuming the temperature on 
the left to he Tl = T + l/2Ar, and on the right, 
Tr = T - l/2Ar, with Ar,w„ < T, we can obtain 
both the heat flux J in the system, and the temperature 
profile T„: 



rriN 

T^ = T+lAT{In~In), 
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where the quantities 
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I-n = '^^^ / d(jJUj'^\GnN? 



(10) 



are expressed in terms of the Green's function: 

G{u:) = {K - Mu? - T.l{uj) ~ . 

with M being a diagonal matrix with elements cor- 
responding to the particle masses, and K, the dynamic 
matrix for the system. The function S(a;) specifies the 



Re / I3{x)dx = k{uj) 
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In particular, for the case of linear mass distribution ^ 
we come to: 



fc(w) = Re 



/(/3(0))-/(/3(l)) 
A 



where /(/3) = (sin/3 — /3cos/3). This dispersion relation 
is displayed in Fig. [2] This result correlates with numer- 
ical simulations [13_. 

As seen from Fig. |2] in systems with linearly in- 
creasing masses one can identify two kinds of vibrations, 
more specifically, delocalized phonons with frequencies 
UJ < V^oJo ■ and localized "gradons" with frequencies 

\/2uJo <UJ < V^UJi, Wo,! = a/2X/Mo,1. 

In a general case, there exist also phonons with fre- 
quencies UJ < \plujmin and gradons with frequencies 
w > \/2uJmin, with uJmin being determined by the maxi- 
mum mass in the system. Incidentally, for gradons /? be- 
comes imaginary, with the determinant starting to grow 
exponentially. As a result, frequencies uj > \/2uj„im do 



not contribute to the integrals in Eqs. (10 I 



Standard methods of averaging [2 permit reducing the 
integrands in Eq. (llOl) in the N oo limit to the form 
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FIG. 2. Dispersion relation of the vibration energy in a har- 
monic chain with linearly increasing masses ([5| plotted vs. 
the wave vector for different relations between the boundary 
masses cjo.i ~ \/2/ir/Mo,i. For energies uj < \/2ujq, the sys- 
tem operates with delocalized phonons, and for \/2cjo < u < 
\/2i^i, with localized gradons. 



mal insulator, J ~ N~'^/'^, and for s = the flux does 
not depend on system size. In the intermediate case of 
s — 2/3, the system conforms to Fourier's law J ~ . 
Note that the dependence of heat flux on the size of a 
harmonic system is specifled by boundary conditions [12J, 
and it is possible to select them such that Fourier's law 
in the system will be obeyed. In the s = 2/3 case we 
selected appropriately the system itself while leaving un- 
changed standard boundary conditions. Thus a system 
with linearly increasing masses can be adjusted such that 
it will conform to Fourier's law. 

Let us turn now to the temperature proflle. First, as 
seen from relations ( 14 1 , if the masses are not constant 



within the chain, the temperature at the points at which 
M'{a) = should rise strongly. 

Second, it can be shown that in the case of constant 
masses thermal proflle within the system T(a) = T is 
constant, a point well enough known. In the reverse ap- 
proximation, in the case of linearly increasing masses, 
when Ml ^ Mq, we come to: 
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Ta^T AT arcsin (Va) ■ 



where n is used in the form n = aN. To preclude misun- 
derstanding, we note that these relations do not go one 
into another at a = 1 and a — 0, because TV — cx) is a 
non-uniform limit. 

Putting relations (14) into we come to a general 
answer for the heat flux and temperature proflle in a 
system with an arbitrary mass distribution. While these 
integrals cannot be performed in a general case, they per- 
mit a number of conclusions. 

To begin with, it turns out that a flow in a system 
depends only on the boundary masses and the maxi- 
mum mass in the chain (which specifles the minimum 
frequency). We readily see that if the masses are con- 
stant, we come to the standard answer for the heat flow 
[3]. In the reverse approximation, when Mi ^ Mq, we 
obtain: 

^ 2V2AT f4>(0)V f r- , /K ~.\ 

J = 5- — — V2ujminl - arctan(V2a;,„i„7 

Up to this point, we assumed the boundary masses not 
to depend on N. We can include this dependence in the 
case of linearly increasing masses. Setting Mi ~ MqN'^ , 
the heat flux will acquire the form: 
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It is easy to verify the correctness of the assumptions we 
have made here, because we know the exact solution for 
the linear case. 

As follows from the expression for the flux J, for a 
constant gradient s = 1 the system behaves as a ther- 



We see immediately that the temperature in the system is 
determined by the more massive end. Now at the center 
of the system under consideration a temperature gradient 
— AT/tt will appear. 

In conclusion, we have presented the solution to 
systems of linear equations ^ describing systems of 
functional-gradient materials. In particular, we analyzed 
the problem of a nonequilibrium steady state of a har- 
monic variable-mass system connected with heat reser- 
voirs which are maintained at different temperatures. 

We have shown that in the particular case of linearly 
increasing masses the heat flux depends on system size 
J ^ with the exponent in this relation being de- 

termined by the law governing the increase of the bound- 
ary mass M ~ iV^. This result flnds ready explanation 
when we turn to Fig. |2] The systems under consideration 
maintain vibrations of two types, delocalized phonons 
which transport heat ballistically and localized "gradons" 
. By properly varying the boundary mass as a function 
of N, we modulate in this way the number of phonons, 
to arrive flnally at the dependence of the heat flux J on 
the system size N. A similar effect accounts for the ap- 
pearance of a thermal gradient in a system with linearly 
increasing masses. We have demonstrated in the partic- 
ular example of a system with linearly increasing masses 
that a harmonic system can both sustain formation of a 
temperature gradient and conform to Fourier's law. 
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tance he has been rendering throughout my work on the 
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